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Thermodynamic quantities in symmetric nuclear matter are calculated 
in a self-consistent T— matrix approximation at zero and finite tempera- 
tures. The internal energy is calculated from the Galitskii-Koltun's sum 
rule and from the summation of the diagrams for the interaction energy. 
The pressure and the entropy at finite temperature are obtained from the 
generating functional form of the thermodynamic potential. 
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1. Introduction 

Nuclear matter is studied as a many-body system of interacting fermions. 
Its equation of state is of particular interest because of its applications to 
heavy-ion collisions, supernovae collapses and neutron stars structure. The 
thermodynamic properties of such systems can be investigated starting from 
the nucleon-nucleon potential; strong interactions between nucleons induce 
short range correlations in the dense nuclear matter. 

Our approach is based on finite-temperature Green's functions. We 
adopt the T— matrix approximation PP, which consists in expressing the 
self-energy as a sum of ladder diagrams. Such a self-consistent scheme 
belongs to a class of approximations derivable from a suitably chosen gener- 
ating functional which automatically fulfill thermodynamic relations, 
including the Hugenholz-Van Hove and Luttinger identities [3J E] • 
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For what concerns the binding energy at zero temperature the thermody- 
namically consistent T— matrix approximation [31 H3 yields results similar to 
other approaches, such as extensive variational and Brueckner-Hartree-Fock 
calculations |Hj • It is hazardous to compute the pressure and the entropy 
simply by integrating thermodynamic relations because of uncontrollable 
numerical errors that propagate at small temperatures. We present instead 
a method which allows to reliably calculate the pressure (and consequently 
the entropy) at sufficiently high temperatures from the diagrammatic ex- 
pansion of the thermodynamic potential [Sj. 

The starting point of these studies is a free N-N potential, a parame- 
terization of the interaction derived from nucleon-nucleon scattering experi- 
ments. In the present paper we restrict ourselves to the two-body CD-Bonn 
potential and will consider three-body interactions in a future work. 



2. T— matrix approximation 

In the in-medium T— matrix approach the two-particle propagator takes 
the following form (for simplicity we skip spin and isospin indices, as well 
as energy and momentum dependence) 

Q 2 = G" c + G" c T G% c + exchange terms . (1) 

G~2 C is the non correlated two-particle Green's function, just constructed as 
a product of two dressed single-particle propagators. We remark that the 
use of the dressed propagator implies the presence of a nontrivial dispersive 
self-energy which leads to a broad spectral function ^U]. The in-medium 
two-particle scattering matrix T, which appears in the above expression, is 
defined as 

T = V + V G% c T (2) 

where V is the interaction potential. All the ingredients are then calculated 
iteratively: the scattering matrix, the single-particle self-energy, which is 
expressed in terms of the scattering matrix, and the single-particle propa- 
gator, by means of the Dyson equation. 

The T— matrix approximation scheme can be as well derived from a 
generating functional V] which is a set of two-particle irreducible dia- 
grams [21 obtained by closing two-particle ladder diagrams with a combina- 
torial factor l/2n, where n is the number of interaction lines in the diagram. 
The functional <3? depends on the dressed propagators (lines in Fig. 0) and 
on the two-particle potential (wavy lines in Fig. 0). The self-energy is then 

5$ 

expressed as a functional derivative S = ■—=. 
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Fig. 1. Two-particle irreducible functional for the T— matrix approximation. The 
sum runs over the number of interaction lines n in the two-particle ladders. 



3. Properties of the correlated system 

3.1. Internal energy 

The (total) internal energy per particle can be calculated as the expec 
tation value of the Hamiltonian of the system H = + H po t 
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Here p is the density and V the volume of the system. When we evaluate 
the right hand side of (jSJ) in momentum space we see that the kinetic part 
is 

where A(jp,u) is the spectral function and f{oj) the Fermi-Dirac distribu- 
tion. For the potential term, using the approximation (^Q) for the two- 
particle Green's function, we get 

V f d 3 P d 3 k dn 
{Hpot) ~ 2 J (2vr)3 ( 2 ,)3^ 6 ^ 

xIm{((k|T i? (P,^)|k)-(k|T iJ (P,0)|-k}) GT R (P,K^)} , (5) 

where b(£l) is the Bose-Einstein distribution. The dependence on the total 
energy £1 and momentum P of the pair and on the exchanged momenta k 
is now made explicit. 

It is also possible to determine the energy in an alternative and simpler 
way. The Galitskii-Koltun's sum rule [TT1IT21IT3] 
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for conserving approximations is equivalent to the direct calculation ([5]). 
However this is valid when only two-body interactions are considered. In 
the presence of many-body forces the sum rule cannot be employed and 
the internal energy of the correlated system has to be computed from the 
expectation value of the Hamiltonian jSJ). 

We calculate the internal energy per particle with the use of the two 
methods: from the Galitskii-Koltun's sum rule © and from diagram sum- 
mation, i.e. from eq. © together with Q) and ©. We restrict ourselves 
to the empirical saturation density (p = 0.16 fm~ 3 ) and consider different 
temperatures up to 20 MeV. The calculation are performed using a numer- 
ical procedure where the energy range is limited to an interval [— uj c ,u; c ]. 
We have tested several values of the energy scale lo c between 2 GeV and 
10 GeV, and we found that the Galitskii-Koltun's sum rule expression for 
the energy shows some dependence on this value. On the other hand, the 
results obtained from the direct estimation of the interaction and kinetic 
energies are stable and independent on the energy range taken (up to an 
inaccuracy due to numerical discretization of about 0.5 MeV). The cutoff 
dependence is shown in Fig. |2]for two different temperatures. The result 
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Fig. 2. Internal energy per particle as a function of the cutoff at T = and 
T = 10 MeV. The dashed lines represent the energy calculated from the sum rule 
©, the solid lines the energy from the summation of diagrams. 
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from the diagram summation can be compared to Galitskii-Koltun's sum 
rule result extrapolated to infinite energy range (Table^). This extrapolated 
value will be used in the following. 

3.2. Pressure 

The pressure is related to the thermodynamical potential through 

n(T,fi,V) = -PV. (7) 

It can be shown that 

n = -Tr{ln[G- x ]} - Tr{£G} + $ . (8) 

The first two terms involve an integration over energy and momenta of 
single-particle propagators and self-energies. The calculation of the func- 
tional $ requires instead a summation of a set of diagrams, evaluated in 
the following way. One notices that within the T— matrix approximation 
the diagrammatic expansions for the interaction energy and the functional 
differ by the factor 1/n where n is the number of interaction lines in the 
diagram . Hence the functional can be obtained from the formula for the 
interaction energy 

rl AX 

<*>= / -T- < H pot (\V,G x=1 ) > . (9) 

JO A 

In the above formula the interaction potential is multiplied by the factor A 
but the propagator G is the dressed nucleon propagator corresponding to 
the system with the full strength of interactions (A = 1). 
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Egk/N 


Ediag/N 


Ptot/P 


Pbf/p 





-15.80 


-16.63 


-7.69 


-5.58 


2 


-15.15 


-16.29 


-5.86 


-5.4 


5 


-14.40 


-15.24 


-5.48 


-5.2 


10 


-11.15 


-11.72 


-2.66 


-3.35 


20 


-1.29 


-1.21 


6.49 


4.74 



Table 1. Results (at p = po) of the internal energy and the pressure at zero and 
finite temperature (all values in MeVs). The second and third colum represent the 
internal energy per particle obtained from the Galitskii-Koltun's sum rule JHJ and 
from the diagram summation J3J . The other columns contain our estimate for the 
pressure (from eq. Q) and the results of Baldo and Ferreira |15j. 

The pressure in hot nuclear matter has been obtained using two-body 
Argonne V14 interaction in the Bloch-De Dominicis approach ^Sj. The 
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results are qualitatively similar, with a negative value of the pressure at 
T = and p ~ 0.16 fm~ 3 . This shows the need to include three-body forces 
for a reliable description of the thermodynamics of the nuclear matter. 

3.3. Entropy 

We compute the entropy through the thermodynamic relation 



S_ _ 1 
N ~ T 



E , Ptot 

U 

N p p . 



(10) 



The results are plotted in Fig. EJ The error in the calculation of TS at each 
temperature can be estimated by comparing the results obtained with the 
two expressions for the internal energy © and © ■ The difference is of the 
order of 1 MeV, also at zero temperature we find ITS'! ~ lMeV ^ 0. The 
entropy can be estimated reliably only for T > 5 MeV, with the uncertainty 
shown as the hatched band in the figure. 

We compare these results with other methods of calculating the entropy: 

1. The dynamical quasiparticle formula |16j 
Sdq 1 f d 3 p duj 



N pj (2?r) 3 2tt 



SRe£*(p,u;) 
A(p,u) 1 - - 



+ r(p ' w) 



in 



where 



a(u) = -f(u) ln[/H] - [1 - f(u)] ln[l - /(c)] . (12) 



It has been shown that this one-body formula gives results close to 
the complete expression for the entropy |17j . 

The entropy calculated as for a free Fermi gas but using the effec- 
tive mass m* (determined at each temperature by {duj p / dp 2 ) p=pp = 
l/2m*) instead of the rest mass m 



,2 



Sfree* _ 1 f d 3 p ( p_ 
N ~ pj (2vr) 3 ° I 2m* 



(13) 



The entropy estimated by means of (Jll[) and (|13|) is shown in Fig. |3]with the 
dashed and solid line respectively. Remarkably, we find that the expression 
for the entropy of the free Fermi gas is very similar to the result of the 
full calculation for the interacting system, if the change of the effective 
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Fig. 3. Entropy per baryon as a function of the temperature. The hatched band 
denotes our estimate of the entropy. The solid line denotes the result for the free 
Fermi gas with the in-medium effective mass and the dashed line represents the 
result of [T7j (Eq.[ITJ). 



mass in the system is taken into account. The last observation simplifies 
significantly the modeling of the evolution of protoneutron stars |18| . since 
relations between the entropy per baryon and the temperature derived for 
a fermion gas can be used in hot nuclear matter. As observed in ^7] the 
quasiparticle expression (fTT|) for the entropy ^H] follows closely the full 
result. 

4. Summary 

We investigate the properties of correlated nuclear matter up to T = 
20 MeV. We calculate the internal energy with different methods, through 
the Galitskii-Koltun's sum rule and by summing the diagrams which con- 
tribute to the expectation value of the interaction energy. The two methods 
give similar results at all temperatures up to a difference of about 1 MeV 
which can be attributed to numerical inaccuracies. The pressure is derived 
from the thermodynamic potential which is given by the generating func- 
tional The diagrams contributing to $ are calculated with an integration 
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of the interaction energy over an artificial parameter A which multiplies the 
interaction lines, while keeping the propagators dressed as in the fully cor- 
related system. We finally estimate the entropy and compare it to two 
expressions derived in the context of a quasiparticle gas, whose results turn 
out to be close to the full calculation. 

Only a study which embodies three-nucleon interactions can be realisti- 
cally compared to experimental data on the nuclear equation of state. Here 
all calculations are carried out at the empirical saturation density, but the 
present scheme can be applied to other densities and it is suitable for the 
inclusion of three-body forces or the study of asymmetric nuclear matter, 
which will be addressed in a future work. 
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